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ABSTRACT 

In  this  paper,  we  study  the  performance  of  the  plasma  opening  switch  (POS)  and  electron-beam  diode  load 
on  the  DM1  pulsed  power  generator  using  the  transmission-line  circuit  code  BERTHA,  the 
electron/photon-transport  code  CYLTRAN,  and  the  2-dimensional  particle-in-cell  (PIC)  code  MAGIC. 

The  results  suggest  that  of  the  52  kJ  required  to  be  delivered  to  the  load  to  account  for  the  observed 
radiation,  up  to  75%  is  associated  with  current  not  flowing  in  the  cathode.  This  non-cathode  current 
cannot  be  explained  by  electron  vacuum  flow  only,  implying  the  presence  of  plasma.  Further,  the  inferred 
experimental  diode  impedance,  1 1  fi,  is  much  higher  than  optimum.  Assuming  no  change  in  the  DM1 
conduction  current/time  and  POS  performance,  an  optimum  diode  impedance  of  2  Q,  and  a  POS-to-load 
inductance  reduced  by  75%,  the  analysis  suggests  that  the  load  energy  on  DM1  can  range  from  54  kJ  to 
94  kJ,  depending  on  how  much  non-cathode  current  reaches  the  load.  The  analysis  also  suggests  that  to 
obtain  agreement  between  features  of  the  measurements  and  calculations,  the  DM1  flow  impedance  (gap) 
must  decay  after  opening  and  ions  must  be  present  in  the  vacuum  transmission  line  between  the  POS  and 
load.  Understanding  of  the  coupling  between  the  POS  and  e-beam  load  remains  incomplete. 

INTRODUCTION 

DECADE1  is  a  proposed  9-MJ,  modular  pulsed  power  generator.  One  such  module,  DECADE  Module  1 
(DM1),2  is  now  in  operation  with  an  electron-beam  (e-beam)  diode  load.  Each  module  uses  a  coaxial  plasma 
opening  switch3  (POS)  for  power  conditioning.  In  this  paper,  we  analyze  the  performance  of  DM1  during  a 
selected  experimental  series.  The  purpose  of  this  work  is  to  try  to  improve  the  understanding  of  POS/load  coupling 
and  to  suggest  ways  to  optimize  the  DM1  single  shot  performance.  Three  different  tools  are  used  in  this  study. 

The  transmission  line  code  BERTHA4  is  used  for  circuit  modeling  that  links  the  pulsed  power  to  the  POS  and, 
using  physics-based  circuit-element  models  for  the  POS  and  e-beam  diode  load,  predict  the  energy  coupled  to  the 
load  (for  the  experiments  considered  here,  there  are  no  load-current  and  only  indirect  load-voltage  measurements). 
The  coupled  electron/photon  Monte  Carlo  transport  code  CYLTRAN5  is  used  to  link  the  measured  radiation 
associated  with  the  e-beam  load  to  the  energy  delivered  to  the  load.  The  2-D  PIC  code  MAGIC6  is  used  to  help 
understand  the  electron  and  ion  flow  between  the  POS  and  load  that  accompanies  the  POS  opening. 

A  simplified  representation  of  the  DM1  circuit  and  an  illustration  of  the  front-end  geometry  are  given  in 
Fig.  1.  An  in-depth  description  of  the  DM1  circuit  can  be  found  elsewhere.1,7  The  front-end  geometry  shown  in 
Fig.  1  was  used  in  the  selected  DM1  experiments  and,  unless  otherwise  noted,  in  all  the  calculations  that  follow. 
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Fig.  1  Representation  of  the  DM1  circuit  and  sketch  of  the  front-end  geometry  used  in  the  experiments. 
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RESULTS  FROM  CIRCUIT  MODELING 

The  48-element  circuit  model  used  in  BERTHA  for  these  circuit  calculations  has  been  successfully 
benchmarked  against  other  models  and  successfully  (within  10%)  compared  with  measured  DM1  electrical 
waveforms  up  to  the  POS.7  Short  descriptions  of  the  circuit-element  models  for  the  POS  and  e-beam  load  as  well 
as  an  example  of  the  circuit  modeling  results  follow.  Details  regarding  the  circuit-element  models  and 
benchmarking  work,  and  complete  circuit  modeling  results,  can  be  found  elsewhere.7  We  have  assumed  the  full 
vacuum  inductance  for  the  region  between  the  POS  and  load.  This  assumption  is  discussed  later.  The  Marx 
charging  voltage  in  the  circuit  calculations  had  to  be  reduced  to  90%  of  the  actual  charging  voltage  in  all  of  this 
work  to  obtain  agreement  with  measurements.7 

The  circuit-element  model  for  the  POS  is  based  on  the  concept  of  flow  impedance,8  where  the  resistance  of  the 
POS  circuit  element,  Rs  is  related  to  the  flow  impedance,  Zf,  by7 


Rs  -  Z, 
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Here,  IAu  and  Icd  are  the  anode  current  on  the  generator  side  (upstream)  and  cathode  current  on  the  load  side 
(downstream)  of  the  POS,  respectively;  Rs  is  the  resistance  of  the  POS  (Rs  =  Vs/Is,  where  Vs  and  Is  =  Iau  -  Icd  are 
the  POS  voltage  and  current,  respectively),  and  Zf  s  Vs/(I2Au  -  I2cd)1/2  In  Eq.  (1)  we  have  also  assumed  that  all  the 
vacuum  electron  flow  between  the  POS  and  load  is  lost  and  counted  as  switch  loss.  No  vacuum  electron  flow  is 
counted  as  current  available  for  driving  the  e-beam  load.  The  load  current  is  Icd-  Because  vacuum  flow  can 
contribute  to  the  radiation  produced  in  the  experiment,  we  expect  this  POS  circuit  model  to  underestimate  the  load 
current.  The  biggest  advantage  of  this  model  is  that  it  allows  one  to  differentiate  between  cathode  and  anode 
current  reaching  the  load.  However,  Zf  still  must  be  prescribed  as  a  function  of  time.  The  physical  model  for  the 
POS  used  in  this  work  assumes  that  the  vacuum  electrons  are  always  critically  insulated  (saturated  flow)  in  the 
time  dependent  POS  gap,  Dcrjt-3'7  If  we  further  assume  that  the  POS  gap  forms  at  the  POS  cathode  radius,  Rc,  then 
one  has  that  Zf «  30  Dmt/Rc.9 

The  circuit-element  model7  for  the  e-beam  load  is  divided  into  three  phases:  1)  an  early  turn-on  phase,  2)  a 
Child-Langmuir  (C-L)  flow  phase,  and  3),  a  critical-current,  pinched-beam  phase.  The  model  uses  as  input:  the 
diode  cathode  outer  radius,  the  cathode  inner  radius,  the  radius  of  curvature  of  the  cathode  edges,  and  the  initial 
anode-cathode  (AK)  gap,  Dd.  The  parameters  for  the  load  model  are:  the  electric  field  turn-on  threshold,  the  time 
duration  to  reach  full  C-L  current,  and  the  gap  closure  speed,  v.  In  the  first  phase,  the  current  begins  to  flow  in  the 
diode  when  the  mean  electric  field  reaches  a  preset  threshold.  Once  the  field  threshold  is  exceeded,  current  rises 
linearly  from  zero  to  the  full  C-L  value  in  the  assumed  diode  “turn-on”  time  (during  which  plasma  expands  to 
cover  the  entire  cathode  area).  In  the  second  phase,  the  full  C-L  current  flows  and  the  diode  AK  gap  closes  from 
electrode  plasma  motion  at  a  constant  speed  (closure  is  initiated  at  the  start  of  “turn-on”).  When  the  C-L  current 
exceeds  the  critical  current,  the  third  phase  begins  which  limits  the  diode  current  to  the  critical  current  for  the 
remainder  of  the  simulation.  The  critical  current  is  expressed  as 

lour  =*500  -\)m  (A).  (2) 

where  Fcr  is  a  factor  that  is  weakly  dependent  on  voltage  and  is  approximately  1.6  for  voltages  between  1  and 
3  MV,  and  y  is  the  ratio  of  electron  energy  to  electron  rest  energy.  The  values  of  the  diode  parameters  that  best  fit 
the  DM1  data,  and  used  here,  are  10  kV/cm  for  the  electric  field  turn-on  threshold  and  0  ns  for  the  duration  to  full 
C-L  current.  Typical  values  of  the  electric  field  turn-on  threshold  and  duration  to  full  C-L  current  for  similar 
diodes  on  the  conventional  water-line  generators  are  200-300  kV/cm  and  10-20  ns,  respectively.  Long-conduction¬ 
time  POS  data  from  Hawk,  where  POS  plasma  may  be  in  the  e-beam  load  region  at  the  time  of  opening,  exhibit  the 
same  values  as  used  here  for  DM1. 10  For  all  of  this  work,  the  closure  speed  was  chosen  to  be  3.5  cm/ps,  typical  of 
similar  diodes  on  water-line  machines. 

For  the  time  dependence  of  Zf  in  Eq.(l),  we  assumed  a  linear  rise  followed  by  an  exponential  decay.7  This 
decay  is  needed  to  reproduce  the  measured  upstream  waveforms.  We  obtained  a  good  fit  to  DM1  data  when  Zf 
linearly  increased  to  1.38  Q  (maximum  POS  gap  of  about  2.3  mm)  in  28  ns,  then  decreased  exponentially 
(corresponding,  e.g.,  to  POS  gap  closure)  with  a  1/e  time  of  55  ns.  Comparison  of  this  model  with  the  measured 
upstream  waveforms  from  a  DM1  shot  is  given  in  Fig.  2.  Also  shown  in  Fig.  2  is  the  calculated  IAu  with  a  short 
circuit  at  the  POS.  It  appears  that  the  present  DM1  operation  utilizes  only  about  78%  of  this  available  current  or 
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about  60%  of  the  available  stored  energy  if  the  POS  conducted  the  entire  current  pulse.  Displayed  in  Fig.  3  are  the 
calculated  load-voltage,  VL,  and  load-current,  IL,  waveforms  along  with  measured  and  calculated  (using  IV  and 
IV2  8  as  approximations)  x-ray  signals.  We  do  not  have  measured  electrical  load  waveforms  to  compare  with  these 
calculations:  however,  the  peak  load  voltage  of  1.7  MV  (inferred  from  a  multi-channel,  filtered  PIN  detector 
array11)  agrees  with  the  calculation.  The  AK  gap  of  25.4  mm  used  in  the  experiment  corresponds  to  a  load 
impedance  of  ZL  =  1 1  Q.  The  calculated  peak  IL  is  150  kA  and  the  load  energy  reaches  12.4  kJ  when  the  measured 
x-ray  pulse  reaches  zero. 

RESULTS  FROM  CYLTRAN 
MODELING 

A  foil  description  of  the  e-beam  diode 
geometry  depicted  in  Fig.  1  and  used  in 
the  CYLTRAN  simulations,  as  well  as 
detailed  results  of  the  radiation 
calculations,  can  be  found  elsewhere.12 
Briefly,  the  12.7-cm  radius  Ta  used  as  an 
anode  in  the  e-beam  diode  was  0.003 -in 
thick  and  the  vacuum  window  was  a  3/16- 
in  thick  sheet  of  polyethylene  followed  by 
a  3/4-in  slab  of  polyethylene.  The  e-beam 
cathode  radius  of  6.35  cm  was  assumed  to 
be  the  e-beam  radius.  The  incident  e- 
beam  angular  distribution  was  forward 
11  12  13  14  15  16  directed.  The  calculated  dose-area- 

Time  fos)  product  (DAP)  in  the  near  field  is 

Fig.  2  Comparison  of  measured  and  calculated  upstream  waveforms,  relatively  insensitive  to  this  distribution. 

Figure  4  contains  a  plot  of  the 
CYLTRAN-calculated  DAP  in  20  mils  of 
Si  as  a  function  of  diode  energy  for  a 
variety  of  cases.  The  horizontal  line  at 
7.6  Mrad  (Si)-cm2  is  the  average 
measured  DAP  for  the  chosen  DM1  data 
set.  The  standard  deviation,  maximum, 
and  minimum  DAP  values  are  also 
shown.  The  open  symbols  are  CYLTRAN 
calculations  for  triangular-voltage,  flat- 
current  load  waveforms.  For  the  circles, 
the  voltage  is  fixed  at  1.75  MV  and  the 
current  is  varied  to  generate  the  linear 
scaling  of  calculated  DAP  as  a  function  of 
load  energy.  For  the  two  types  of 
triangles,  the  current  is  held  fixed  at 
either  0.6  MA  or  1.2  MA  and  the  voltage 
is  varied  between  1  and  2  MV.  Note  that 
Fig.  3  Calculated  load  parameters  and  measured  x-ray  signal.  for  these  geometries  the  DAP  is  an 

approximately  linear  function  of  load 

energy.  This  result  means  that  electrical  energy  delivered  to  the  diode  is  a  good  parameter  for  optimizing  radiation 
production.  The  filled  squares  are  calculated  DAPs  for  various  BERTHA-generated  load  waveforms.7  The  point 
labeled  “Dcm(t)  model”  corresponds  to  prescribing  Zf  in  Eq.  (1)  by  calculating  Dcnt(t)  from  a  form  of  Eq.  (2)  using 
the  measured  upstream  waveforms  and  setting  Zf  =  30  Deft/Rc-7  The  point  marked  “100  ns”  is  the  result  when 
only  the  energy  in  the  first  100  ns  of  the  calculated  power  pulse  is  used.  More  energy  is  delivered  to  the  load  when 
the  AK  gap  and  switch-to-load  inductance,  LSl,  are  reduced  because  the  present  system  is  not  optimized  from  a 
circuit  point  of  view  (see  next  section).  The  CYLTRAN  results  for  the  BERTHA  waveforms  fall  on  the  same  line 
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as  the  idealized  waveforms,  indicating  that  the  DAP  in  20  mils  of  Si  for  this  configuration  is  relatively  insensitive 
to  the  shape  of  the  voltage  and  current  waveforms. 


The  energy  at  which  the 
calculated  DAP  vs.  load  energy  curve 
intersects  the  average  measured  DAP 
is  the  energy  we  infer  is  producing 
the  observed  DAP  on  DM1  -  about 
52  kJ.  The  BERTHA  calculation 
underestimates  the  energy  delivered 
to  the  load  by  approximately  40  kJ. 
The  circuit-element  model  for  the 
POS  assumes  only  current  flowing  in 
the  metal  cathode,  150  kA  peak, 
reaches  the  load.  Thus  we  conclude 
that  in  the  experiment,  up  to  40  kJ  of 
energy  associated  with  non-cathode 
current,  525  kA  (assuming  the 
experimental  load  power  pulse 
waveform  is  comparable  with  the 
calculated  one),  reaches  the  anode 
and  makes  radiation.  Electron 

Fig.  4  Dose-area-product  as  a  function  of  diode  energy.  vacuum  flow  310116  cannot  exPlain 

this  result  (see  below).  Note  that 

comparisons  of  measured1 1  and  calculated  radial  dose  profiles  suggest  that  many  electrons  hit  the  anode  at  a  radius 
well  outside  the  e-beam  cathode  radius,  supporting  our  conclusion  that  non-cathode  current  is  present  in  the 
experiment. 


RESULTS  FROM  VARYING  CIRCUIT  PARAMETERS 

Several  BERTHA  runs  were  made  varying  Dd  and  LSl  from  25.4  mm  and  63  nH  to  5  mm  and  15.8  nH, 
respectively  (the  former  values  being  the  ones  used  in  the  experiments).7  These  calculations  show  that  the  energy 
delivered  to  the  load  has  a  broad  maximum  of  37  kJ  when  Dd  is  about  6.3  mm,  with  an  associated  ZL  of  about  2  Q. 
This  optimization  was  expected  based  on  the  critical  gap  picture  of  POS  opening,  where  the  optimum  load 
impedance  for  maximum  energy  delivery  to  the  load  defines  the  boundary  between  switch-limited  and  load-limited 
operation.3 

If  Lsl  is  reduced,  the  load  energy  is  also  increased,  but  only  if  the  load  impedance  is  near  the  optimum 
impedance.  If  ZL  is  large  compared  to  the  optimum  ZL,  the  beneficial  effect  of  reducing  LSl  is  negated.  For 
example,  with  Dd  =  25.4  mm,  there  is  no  change  in  the  energy  delivered  to  the  load  as  Lsl  is  reduced.  This  is 
because  ZL  is  so  large  (1 1  Q)  it  dominates  the  inductive  loading  (approximated  by  Lsi/topen ,  about  2  Q  for  LSl  of 
63  nH  and  t^n  of  28  ns).  ZL  of  about  2  Q  is  comparable  to  the  maximum  inductive  loading  and  the  system  can 
deliver  more  energy  to  the  load  as  Lsl  is  reduced.  When  Dd  =  6.3  mm  and  LSl  =  15.8  nH  (25%  of  the  original 
vacuum  inductance),  we  calculate  that  54  kJ  is  delivered  to  the  load.  If  we  conjecture  that  the  40  kJ  of  load  energy 
associated  with  the  non-cathode  current  is  still  delivered  to  the  load,  then  up  to  94  kJ  (upper  bound)  might  be 
delivered  to  the  load  with  this  front-end  configuration.  The  real  system  response  must  be  determined 
experimentally. 


RESULTS  FROM  2-D  PIC  MODELING 

To  gain  insight  into  the  details  of  the  POS-load  coupling,  we  carried  out  preliminary  2-D  PIC  calculations  of 
the  behavior  of  the  ion  and  electron  flows  using  MAGIC.  The  DM1  front-end  geometry  (see  Fig.  1)  was  used  with 
the  POS  modeled  by  independently  moving  ions  out  of  a  prescribed  gap  in  a  prescribed  time,  as  described 
previously.9  For  the  results  shown  here,  the  final  gap  was  1.8  mm  and  the  time  for  ions  to  clear  the  gap  was  30  ns. 
This  model,  together  with  a  simplified  driving  circuit,  gave  a  good  match  to  the  upstream  DM1  waveforms.  POS 
plasma  was  modeled  by  allowing  ion  emission  along  both  sides  of  the  POS.  The  load  for  the  simulation  was  an  e- 
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beam  diode  with  no  ions,  an  electric  field  threshold  for  electron  emission  of  100  kV/cm,  and  a  25-mm  AK  gap. 

The  full  DM1  current  flowed  in  the  circuit  upstream  of  POS  at  time  t=0  (initiation  of  opening). 

Measurements  on  Hawk  with  downstream  current  monitors  in  the  anode  and  cathode  suggest  that  the  current 
propagates  to  the  load  nearer  the  anode  first  with  a  speed  of  the  order  of  1000  cm/ps,  followed  by  current  flowing 
in  the  cathode.13  Recent  results  from  DM1  show  similar  behavior.  On  Hawk,  the  time  delay  between  the  current 
measured  near  the  load  and  at  the  anode  and  the  initiation  of  cathode  current  flowing,  as  well  as  the  ratio  of  anode 
to  cathode  current,  decrease  with  decreasing  ZL.  Several  possible  mechanisms  have  been  proposed  to  explain  the 
process  of  current  transfer  to  the  load  including:  vacuum  electron  flow,  electron  migration  along  the  downstream 
anode  via  plasma  production  at  the  anode,14  and  various  scenarios  involving  electron-propagation/current-transport 
through  or  by  a  downstream  plasma.  In  the  simulations  described  here,  there  is  no  attempt  to  include  plasma  that 
may  be  in  the  POS-to-load  region  as  a  result  of  the  POS  plasma  injection  and/or  the  conduction  process.  Also, 
there  was  no  self-consistent  algorithm  for  producing  ions  at  the  anode  from  electron  energy  deposition.  Within  the 
context  of  these  assumptions,  the  POS-to-load  inductance  is  self-consistently  accounted  for  in  the  PIC  simulations. 

In  Fig.  5,  ion  and  electron 

electron  positions  at  t=25  ns  electron  positions  at  t=50  ns  are  plotted  at  ^ 

different  times.  At  t=25  ns,  the 
electrons  fill  the  downstream 
region  and  drag  ions  associated 
with  the  POS  plasma  into  the 
vacuum  line  between  the  POS 
and  load.  At  t=50  ns,  the  ions 
(plasma)  also  fill  the  vacuum 
line.  The  behavior  of  the 
particle  flow  is  distinctly 
different  from  that  reported  in 
previous  simulations,9  where  the 
electrons  hugged  the  anode  of 
the  vacuum  line  between  the 
POS  and  load.  In  this  case,  the 
electron  flow  is  more  uniformly 
distributed  in  radius.  The  cause 
of  this  discrepancy  is  under 

investigation.  The  onset  of  significant  (anode)  current  in  the  diode  is  associated  with  the  arrival  of  plasma  to  the 
diode.  The  speed  with  which  the  plasma  moves  depends  on  the  ratio  of  ion  charge  to  atomic  number,  Z/A,  which 
for  this  simulation  was  assumed  to  be  2/5  (corresponding,  e.g.,  to  an  18%  C++/82%  H+  plasma  -  by  comparison, 
purely  C ",+  ions  would  have  Z/A=2/12).  This  Z/A  was  chosen  to  have  the  anode  current  and  load  power  peak  near 
the  expected  time  after  opening.  The  cathode  current  predicted  by  this  simulation  is  close  to  the  150  kA  predicted 
by  the  circuit  modeling.  Also,  the  sequencing  of  anode  and  cathode  monitors  observed  on  Hawk13  and  DM1  is 
qualitatively  reproduced,  with  the  speed  of  propagation  about  1000  cm/ps.  However,  the  total  current  reaching  the 
load  in  the  simulation  was  only  200  kA,  a  factor  of  3.5  lower  than  what  was  inferred  above  from  the  radiation 
measurements.  With  a  Z/A  of  2/12,  the  total  load  current  is  less  than  150  kA.  These  results  will  depend  on  the 
plasma  Z/A  and  the  existence  of  other  plasma,  either  injected  into  or  generated  in  the  vacuum  line.  For  example, 
in  a  recent  simulation  ion  emission  was  allowed  along  the  anode  (treated  as  a  zero  work  function  ion  source) 
between  the  POS  and  load.  This  resulted  in  a  peak  load  current  of  450  kA,  a  closer  match  to  the  inferred  load 
current,  but  with  the  current  beginning  15  to  20  ns  too  early.  Here,  a  model  for  anode  plasma  production  based  on 
local  energy  deposition  is  needed.  Additional  simulations  along  with  better  measurements  of  plasma  density  are 
required  to  more  fully  understand  the  POS/load  coupling. 
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Fig.  5.  Electron  and  ion  positions  at  two  times  after  opening. 


SUMMARY 

We  have  performed  an  analysis  of  a  selected  set  of  data  from  the  DM1  pulsed  power  generator  using  the 
transmission-line  circuit  code  BERTHA,  the  electron/photon-transport  code  CYLTRAN,  and  the  2 -dimensional 
particle-in-cell-code  MAGIC.  The  results  suggest  that  up  to  75%  of  the  energy  delivered  to  the  diode  is  associated 
with  non-cathode,  radiation-producing  current.  For  the  experimental  series  considered,  the  system  is  not 
optimized  for  energy  delivery  to  the  load.  The  diode  impedance  of  1 1  D.  is  much  higher  than  optimum.  With  the 
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same  DM1  conduction  current/time  and  POS  performance  considered  in  this  paper,  an  optimum  diode  impedance 
of  2  Q,  and  a  POS-to-load  inductance  reduced  by  75%,  this  analysis  suggests  the  load  energy  on  DM1  can  range 
from  54  kJ  to  94  kJ,  depending  on  how  much  non-cathode  current  reaches  the  load.  To  obtain  agreement  between 
features  of  the  measurements  and  calculations,  the  DM1  flow  impedance  (gap)  must  decay  after  opening  and  ions 
must  be  present  in  the  vacuum  transmission  line  between  the  POS  and  load.  While  progress  has  been  made  in 
describing  the  coupling  between  the  POS  and  e-beam  load,  both  experimental  and  theoretical  work  is  required  for  a 
fuller  understanding  of  this  complex  process. 
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